Black hole constraints on the running-mass inflation model 
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The running-mass inflation model, which has strong motivation from particle physics, predicts 
density perturbations whose spectral index is strongly scale-dependent. For a large part of parameter 
space the spectrum rises sharply to short scales. In this paper we compute the production of 
primordial black holes, using both analytic and numerical calculation of the density perturbation 
spectra. Observational constraints from black hole production are shown to exclude a large region 
of otherwise permissible parameter space. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Particle physics models of inflation based on super- 
gravity theories are plagued by the so-called 77-problem 
[pj , which states that the mass-squared of any scalar field, 
including the putative inflaton field, is typically of order 
H 2 (H being the Hubble parameter) which ruins slow-roll 
inflation. An elegant proposal to circumvent this is the 
running-mass model of inflation, introduced by Stewart 
[^],[| , where the flatness of the potential arises because of 
the quantum corrections, which serve to flatten the po- 
tential over a significant region where inflation can then 
take place. 

Because the flatness is brought about by a cancella- 
tion of the intrinsic curvature of the potential against the 
quantum corrections, only a limited portion of the poten- 
tial can support slow-roll inflation, as is necessary to gen- 
erate the approximately flat power spectrum seen by the 
COBE satellite. Well away from COBE scales, one typi- 
cally expects to see dramatic deviations from near scale- 
invariance as the slow-roll regime breaks down. Copeland 
et al. Ql examined the possibility that this breakdown 
of scale-invariance might be detectable through scale- 
dependence of the spectral index of primordial perturba- 
tions, and more recently the model has been extensively 
explored by Covi, Lyth and collaborators |3fj7| in a series 
of papers investigating its viability both from a theoret- 
ical standpoint and in confrontation against large-scale 
structure data. 

In this paper we examine constraints arising from the 
more radical departures from scale-invariance which may 
take place towards the end of inflation. In much of pa- 
rameter space, the spectrum rises sharply on short scales, 
which can give rise to production of primordial black 
holes (PBHs). These are strongly constrained by obser- 
vation and, as we will see, a significant region of otherwise 
viable parameter space is excluded. 
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II. THE RUNNING-MASS MODEL 

Whether or not a potential V(4>) can support slow-roll 
inflation can be judged via the slow-roll parameters S 
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where primes are <p derivatives. When the slow-roll pa- 
rameters are much less than unity, slow-roll inflation can 
proceed and gives rise to perturbations with an approxi- 
mately scale-invariant spectrum. 

Within the context of softly-broken global supersym- 
metry, the false vacuum dominated potential 
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arises naturally [jlj . However it will not in general lead to 
slow-roll inflation, because supergravity corrections lead 
to |/i 2 | = |ry y | ~ 1 in Planck units. In the scenario pro- 
posed by Stewart ^,||, the inflaton has gauge couplings 
to vector or chiral superfields, and one-loop quantum cor- 
rections flatten the potential, corresponding to a running 
of the effective mass with the scalar field value 
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where /i 2 , represents the inflaton mass squared, Aq is 
the mass squared of the gaugino appearing in the loop, 
both evaluated at the Planck scale, and ao is a gauge 
coupling times a group theoretic factor which may be 
positive (in the case of asymptotic freedom) or negative 
(in the opposite case). The functional form of /J, 2 (4>) is 
obtained by solving the relevant renormalization group 
equations. For definiteness, we consider the inflaton po- 
tential |,|,|§ 
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Hubble parameter H and its derivatives, rather than the 
potential. 

The value of rj starts off negative near the maximum 
of the potential and runs through zero until the end of 
slow-roll inflation is reached, which we define to be 



1]{4>end) = 1. 
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At this point the inflaton potential is still dominated by 
the Vq term, so it is assumed that the field must decay 
via some hybrid inflation mechanism when 4> falls below a 
critical value <p c , in order that reheating occurs to restore 
the standard cosmology. 



FIG. 1. A sketch of the potential. The inflaton starts near 
the maximum with r\ v negative. As it rolls towards the origin 
the mass passes through zero and r\ v grows to 1. For the 
purposes of this diagram the bump is greatly exaggerated; in 
reality the potential is extremely flat. 



The term proportional to An vanishes when <f> = Mp and 
the potential reverts to the form of Eq. Far below 
the Planck scale, the desired cancellation occurs between 
the Hq and An pieces allowing slow- roll inflation to occur. 

An attractive feature of this model is that its param- 
eters take on natural values; for ceo > we have Ao and 
//q both positive and of order unity, while for ceo < we 
have Aq negative and the model is subject to the con- 
straint | An | > /Zq + 1, which is applied to ensure that 
the inflaton mass changes sign before reaching the end 
of inflation. Full details of the model can be found in 
Refs. [|§. 

The inflaton field starts off near the maximum of the 
potential^] and rolls towards the origin, corresponding to 
model (i) of Ref. || and as shown in Fig. m. In the 
case of So > 0, the potential has an unphysical pole at 
ln((/>/Mp) = l/a and should not be trusted in this region 
of strong coupling. It is not necessary, though, to evolve 
<f> to such small field values. 

Throughout this calculation we use the Hubble-slow- 
roll parameters, defined as 1 1 1 
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where the fundamental quantity is now taken to be the 



*Note that initial conditions near the maximum of the poten- 
tial are well motivated by the 'topological inflation' idea Q; 
the initial conditions inevitably have the field on different 
sides of the maximum in different regions of space, and hence 
crossing the maximum in the interpolating regions. 



III. COMPUTING THE PERTURBATION 
SPECTRA 

Our main focus is the perturbations near the end of in- 
flation, where the slow-roll approximation will be poor. It 
is therefore imperative that the accuracy of calculations 
is checked numerically. There are two aspects to this; nu- 
merical calculation of the classical scalar field dynamics, 
and numerical calculation of the perturbation equations. 

The numerical calculation of the classical evolution is 
important in determining which part of the potential 
generates the perturbations seen on cosmological scales. 
When the slow-roll approximation begins to break down, 
commonly-used expressions such as that for the number 
of e-foldings N can lose their accuracy and the numeri- 
cal evolution can lead to some corrections to the analytic 
results. In general, N and the wavemode k leaving the 
horizon at that epoch are related by 
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where A^cobe is defined throughout as the number of e- 
foldings before the end of inflation when our present Hub- 
ble radius, in comoving units, equalled the Hubble radius 
during inflation. This expression neglects the variation 
of H, which is valid as long as e v is small even when rj v 
is not. 

To numerically compute the perturbations, we use the 
Mukhanov formalism |ll| , |l2fl as described by Grivell and 
Liddle to which we refer the reader for details. The 
approach involves a numerical solution both of the clas- 
sical homogeneous equations of motion and of the equa- 
tions describing linear perturbations, with the full power 
spectrum being built up mode-by-mode. In order to eval- 
uate the power spectrum, one has to follow the modes 
until they are well outside the horizon, where their ampli- 
tude becomes constant. This becomes problematic once 
one reaches very close to the end of inflation, where this 
asymptotic regime is not reached. In fact, in the case 
of the running-mass model it proves difficult to numeri- 
cally evolve the mode evolution much more than around 
one e-folding beyond the rj — 1 point, since the slow-roll 
parameters rj and £ 2 are growing so rapidly at this point. 
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FIG. 2. An example power spectrum, taking ceo = 0.01, 
Aq — 1.0, fjio — 0.5 and iVcoBE = 45. The numerical cal- 
culation of the perturbation amplitude breaks down towards 
the end of inflation. The end point of the power spectrum 
(fccnd = e NcOBE /i/3000Mpc _1 ) is defined to be where n = 1. 
The position where r\ — 0.5 is shown for comparison. 



be excluded by PBH constraints. The scale-dependence 
of the spectral index is clearly seen. By virtue of the 
N(k) relation of Eq. (Q), a change in A^cobe corresponds 
simply to a translation of the power spectrum to a new 
end point given by k ond = e NcOBE /i/3000 Mpc -1 (though 
it must be renormalized to COBE) . A change in the 
condition for the end of slow-roll inflation also results in 
a translation and renormalization of the power spectrum. 
As we shall see in Section [^L the condition for the end 
of slow-roll inflation, Eq. (g|), is one of two parameters 
that affect the severity of the PBH constraints the most, 
the other being A^cobe- 

In Section IV we consider the cases A^obe = 45, which 



corresponds to instant reheating after the end of slow-roll 
inflation, and A^obe = 25, which may result if the end of 
slow-roll is followed by a bout of fast-roll and/or thermal 
inflation. 



IV. CONSTRAINTS FROM PRIMORDIAL 
BLACK HOLE PRODUCTION 



However, as long as one uses the numerical solution for 
the classical background evolution, it turns out we do not 
need to compute the perturbations numerically; the ex- 
tended slow-roll approximation of Stewart and Lyth 113] 
proves perfectly adequate. This gives the perturbation 
amplitude as 
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where C = -2 + In2 + &~ -0.73 and b is the Euler- 
Mascheroni constant, and gives sufficiently accurate re- 
sults even when rj becomes large. Fig. shows a com- 
parison of numerical simulation with the slow-roll and 
Stewart-Lyth predictions. 

The error in the Stewart-Lyth expression is expected 
to be 0(^ 2 ). It is known to underestimate the perturba- 
tion amplitude when r/ ~ 1 and e«l fl3|| , which makes 
our PBH constraints conservative. A further consequence 
of the smallness of e, 
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is that the gravitational waves from this model are 
strongly suppressed || . 

We use the COBE normalization scheme of Ref. j^J, 
setting £!o = 0.35 and Ha = 0.65. In fact the nor- 
malization of the potential is very nearly independent 
of the cosmological parameters, since the temperature 
anisotropics are, with the exception of the integrated 
Sachs- Wolfe effect, laid down at the redshift of decou- 
pling, long before the cosmological constant is impor- 
tant. In this scheme, normalization occurs at the scale 
k = 7a H = 7/1/3000MPC -1 . 

A typical power spectrum is shown in Fig. ^ taken 
from a region of parameter space that we will show to 



A. The black hole constraint 

The astrophysical details of PBH constraints have been 
studied in detail elsewhere Over a wide range of 

mass scales, the observational constraint on the black 
hole formation rate is that no more than around 10~ 20 
of the mass of the Universe can be channeled into black 
holes. For our purposes it is sufficient to ignore the details 
of the constraints, and simply adopt this level, as the 
black hole production rate is enormously sensitive to the 
amplitude of perturbations. 

In computing the black hole formation rate, the quan- 
tity which is of interest is the matter dispersion a, which 
is defined, in the usual way (see e.g. Ref. ||), as the 
matter distribution smoothed over some length scale i?, 



a 2 {R,t) = 
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We will take W(kR) to be a gaussian filter. The factor 
10/9 appears because we are interested in perturbations 
in the radiation era rather than the usual matter era. 

At the end of slow-roll inflation we have r\ = 1, and 
so the power spectrum is rising steeply with a spectral 
index n — 1 ~ 2rj ~ 2. Approximating b\{k) as a power 
law at the end of slow-roll inflation, 
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and setting aH — 1/R, we can evaluate the dispersion 
Eq. @ at hor izon crossing, 
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FIG. 3. Parameter space constraints for So = 0.01 for two choices of iVcoBE- The solid lines show the spectral index on 
COBE scales. The dashed lines are ow = 0.02, 0.04, 0.08 contours, where ow has been evaluated at the end of slow-roll 
inflation. Parameter space below this region is excluded, and these models will violate the bound on <Th or before the end of 
slow-roll inflation. From Eq. (^), instant reheating requires Vb/Mp > 10 -36 for ./Vcobe = 45, indicated by the dot-dashed 
contour; the parameter space above this contour is excluded. 



where I(n) is a numerical factor of order unity which de- 
pends on the spectral index. The length scale at the end 
of slow-roll inflation also provides the natural scale over 
which to smooth the power spectrum since it is the scale 
on which black holes will predominantly form. Setting 
kcndR = 1 we can evaluate / to be 

r 1 i 

/= / U n+2 ^W 2 Ckk end R)dk = - 7 [(n + 3)/2,l] , (13) 
Jo * 

where 7 [a, x] is the incomplete gamma function. For ex- 
ample, with n — 1 = 3 we have / ~ 0.067, n — 1 = 2 we 
have I ~ 0.080, while for n - 1 = 1 we have I ~ 0.100. 
We note that in the limit of large n, holding <5^(fc cnc j) con- 
stant, the contribution of this spike to Ohor is suppressed 
by the numerical factor /. We can see immediately that 
under the power-law approximation of Eq. (|ll]), the dis- 
persion (Thor only depends weakly on the exact value of 
the spectral index at the end of slow-roll inflation which 
we take to be a nominal and conservative n = 3. 

The main dependence of <7hor is on the perturbation 
amplitude 5n(k cn d), although in our case, and for any 
sharply rising power spectrum, ^(fcond) depends strongly 
on the exact condition for the point where slow-roll in- 
flation ends. 

As shown in Ref . jL7) , the black hole constraint across 
all scales simply amounts to 

tfhor < 0.04 . (14) 

This is sufficient to ensure that no more than 10 -20 of 
the mass density of the Universe is channelled into black 
holes. The constraint on (Thor is expected to be accurate 
to within a factor of 2, which is small compared to the 2- 
3 orders of magnitude that the power spectrum can rise 
between the 77 = 1/2 and rj = 1 points. This uncertainty 
is therefore much less important than the uncertainty of 
the end-point of inflation, which we discuss further below. 



As well as the observational constraints on the model, 
there is a self-consistency constraint which must be satis- 
fied, which is to ensure that the inflationary energy scale, 
once normalized to COBE, is high enough to permit the 
claimed number of e-foldings Acobe- If we conserva- 
tively assume instant reheating after inflation, and that 
the radiation era is not punctuated by episodes of ther- 
mal inflation or temporary matter domination, the upper 
bound on the number of e-foldings that can take place is 

A^cobe < 48 + ln(V 1/4 /10 10 GeV), (15) 

which requires V /M$ > 10~ 36 , 10~ 72 for A CO be = 
45,25 respectively. 

B. Results 

The results for iVcoBE = 45, 25 and <So = 0.01 are 
shown in Fig. We plot three different values of Ohoi-; 
the central one is the best guess at where the constraint 



t If one compares our contours of n on COBE scales with 
those in Ref. ||, differences are apparent especially at large 
values of n. These differences are due to a slightly different 
choice for the end of inflation, and our use of numerical calcu- 
lations rather than an approximate analytic technique. The 
differences should be regarded as indicating the arbitrariness 
in deciding where inflation comes to an end once the form 
Eq. (fj) breaks down away from the slow-roll regime. This 
leads to a different identification of the part of the potential 
corresponding to COBE scales, and the running of n causes 
the contours to slide to a different location. The construc- 
tion of a complete model including a mechanism for ending 
inflation would be needed to remove this uncertainty. 
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FIG. 4. As Fig.^, but for ao = —0.01. The thick line to the right is a bound on the allowed values of the parameters 
{\Aq\ > Ho + 1). The dashed lines, reading from left to right, are the ow = 0.02, 0.04, 0.08 contours, and the region to the right 
of these contours is excluded by PBH constraints. From Eq. (jl^), instant reheating requires Vo/Mp > 10 -36 for -/Vcobe = 45 
indicated by the dot-dashed contour; the parameter space above this contour is excluded. 



lies and the others indicate the uncertainty. We see that 
the PBH constraint actually quite closely follows lines of 
constant n on COBE scales, enabling us to use this to 
summarize the constraint. In the case of Acobe = 45 we 
find that a large amount of otherwise viable parameter 
space is excluded: models with n > 1.1 are ruled out. 
This should be compared with PBH constraints on infla- 
tion models with constant spectral index, for which the 
end result is n > 1.25 are excluded |l6| . It is of course not 
surprising that the constraint on n should be stronger for 
the running-mass model whose spectral index increases 
as a function of wavenumber k. 

For -/Vcobe — 25 the PBH constraints arc less severe, 
models with n > 1.3 being excluded. The simplest expla- 
nation is that the mass runs for fewer e-foldings leading 
to a safer period of slow-roll inflation. In fact it is a com- 
bination of two factors that make the -/Vcobe = 25 case 
safer. Firstly, reducing -/Vcobe has the effect of translat- 
ing the spectral index contours away from the region of 
parameter space previously excluded: for a given spec- 
tral index contour, the values of Aq and /1q must increase 
as -/Vcobe decreases to ensure that the running of the 
mass up to n = 1 is faster. Secondly, when we renor- 
malize the new spectra (for given values of /i§ and Aq) 
to COBE, the overall amplitude at the end of slow-roll 
inflation will be reduced for all models with n > 1.0 as 
compared to the -/Vcobe = 45 case, because the spec- 
tral index n is an increasing function of -V. Therefore, 
since the Chor = 0.04 contour lies in the region n > 1.0 
for A^cobe = 45, the excluded region of parameter space 
shrinks for the A^obe = 25 case. 

Next we look at the results for ao = —0.01, which are 
shown in Fig. [|. For Acobe = 45, the PBH constraints 
are more severe, ruling out models with n > 1.0. This 
result is related to the running strength throughout in- 
flation, given by Eg] (neglecting terms in e) 
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(16) 



Over observable cosmological scales the running strength 
given by — £ 2 is generally greater for negative ao case (for 
a given spectral index contour n) resulting in a larger 
value of n throughout slow-roll inflation. For instance, 
for n — 1 = we have, over observable scales H, 
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Towards the end of slow-roll inflation the running be- 
comes stronger in the positive So case as the field ap- 
proaches our fixed end point, n — 1. For Acobe = 25 
the spectral index contours are once again shifted close 
to the point where the PBH constraints cease to be very 
interesting, with n > 1.3 being excluded. 

Finally we would like to know what effect varying the 
value of ao has on the PBH constraints. Given that the 
Ohor contours are approximately parallel to the spectral 
index contours, ncOBE, we can reduce the dimensional- 
ity of this calculation, and simply assign to each ao some 
critical value of the spectral index (on COBE scales), 
n cr it, above which the model is excluded. We will assume 
the constraint is <7hor < 0.04. The results are shown in 
Fig. |5| and illustrate the trends of the PBH constraints 
when we vary ag. For ao > the PBH constraints be- 
come less restrictive as ao is increased, since the other 
parameters of the model take on smaller values, which 
has the effect of reducing the overall running strength. 
For ao < the PBH constraints become more restrictive 
as | <So| is increased, since the other model parameters 
remain fairly static, and the overall running strength be- 
comes greater. 

Using the so-called linear approximation described in 
Refs. |||7|] , the observational constraints on the running- 
mass model can be expressed in terms of two new pa- 
rameters c and <7, rather than directly in terms of the 



5 



1 20 


1 1 1 1 1 1 1 

5 >0 




a <0 


1.10 




1.00 




0.90 




0.80 





0.01 0.10 
I5 l 



FIG. 5. The effect of varying the coupling constant, ao 
(A^cobe = 45). n C rit is the spectral index on COBE scales 
above which the model is ruled out, assuming the constraint 
is Ohor < 0.04. For ao > the constraint is weakened as ao 
is increased, while for So < the constraint becomes more 
restrictive as |ao| is increased. 

three model parameters /Zq, A and &o- The quantity c 
is related to the coupling strengths involved and a is an 
integration constant related to the endpoint of slow-roll 
inflation. In more general models these two quantities 
are still enough to describe the density perturbation over 
cosmological scales, but not away from these scales where 
the linear approximation breaks down, hence the need 
for a numerical calculation of the PBH constraints. For 
given values of -/Vcobe and n cr ;t, though, the correspond- 
ing constraints on the c-cr plane can be found from Fig. ^ 
using the relation || 

n - 1 = 2ae^ cNcOBE - 2c . (18) 

Before ending, we need to comment on our choice for 
the end of slow-roll inflation given by Eq. (||); as we have 
remarked the constraints can be highly sensitive to this 
and we need to ensure we are being conservative. There is 
no kinematical reason why inflation cannot proceed when 
eCl and i] < 1, although we know that this situation 
can be tolerated for no more than a few e-folds, given 
a COBE-normalized spectrum of perturbations. This is 
just restatement of the 77-problem, and indeed is observed 
in our simulations where inflation always proceeds to the 
7/ = 1 point and beyond. However, as we move into the 
regime where 77 > 1 we find that the running strength 
given by — £ 2 begins to blow up, marking the failure of 
the one-loop approximation. This suggests it is danger- 
ous to try and proceed further along the potential even 
though the numerical simulations show inflation contin- 
uing and the spectrum continuing to rise. Thus, evolving 
the inflaton to </> cn d where r\ = 1 but not beyond ap- 
pears reasonable. Models with a larger spectral index 
on COBE scales will of course violate the bound on <7h or 
before on d is reached and in this way are more strongly 
constrained. 



We are not able to say what happens after the end 
of slow-roll inflation. Eventually inflation is supposed to 
end via the hybrid mechanism when the field passes an 
instability point. However, since the form of the poten- 
tial is breaking down by then we cannot make accurate 
computations in order to check whether there are any 
dangerous perturbations produced during this final era. 
In ignoring such perturbations, we are adopting a con- 
servative approach to the constraints, as our constraints 
from the evolution up to the end of the slow-roll era re- 
main valid whatever might happen subsequently. 

V. CONCLUSIONS 

We have investigated primordial black hole constraints 
within the context of the failure of slow-roll inflation, fo- 
cusing on the well-motivated running-mass model of in- 
flation which features a strong scale-dependence of the 
spectral index 0. Although applying to the amplitude 
of perturbations on very short scales, to a good approxi- 
mation the constraint can be represented as a constraint 
on tt-cobe, the spectral index on the largest observable 
scales. The constraint depends strongly on the number 
of e-foldings -/Vcobe between the production of those per- 
turbations and the end of inflation, and, as with models 
with a constant spectral index, the constraint becomes 
weaker as -/Vcobe is reduced. 

We have shown that a significant region of the pa- 
rameter space of the model, viable under other con- 
straints, is excluded by excess production of black holes. 
This demonstrates the importance of evaluating the den- 
sity perturbation spectrum not just across astrophysical 
scales but also right to the end of inflation. In mod- 
els where the slow-roll approximation holds accurately 
only over a limited range of scales, such as the running- 
mass model, there will be strong deviations from scale- 
invariance towards short scales. In models where the 
deviation takes the form of a strongly blue spectrum, 
excessive black hole production is always likely to be a 
danger. 
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